Buckling Thin Disks and Ribbons with Non-Euclidean Metrics 
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I consider the problem of a thin membrane on which a metric has been prescribed, for example by 
lithographically controlling the local swelling properties of a polymer thin film. While any amount 
of swelling can be accommodated locally, geometry prohibits the existence of a global strain-free 
configuration. To study this geometrical frustration, I introduce a perturbative approach. I compute 
the optimal shape of an annular, thin ribbon as a function of its width. The topological constraint 
of closing the ribbon determines a relationship between the mean curvature and number of wrinkles 
that prevents a complete relaxation of the compression strain induced by swelling and buckles the 
ribbon out of the plane. These results are then applied to thin, buckled disks, where the expansion 
works surprisingly well. I identify a critical radius above which the disk in-plane strain cannot be 
relaxed completely. 

PACS numbers: 46.32.+X, 81.16.Dn,46.70.De 



INTRODUCTION 



In nature, elastic instabilities determine the three- 
dimensional structure of a variety of surfaces, including 
leaves [3, 0, Q , flowers 0| , torn plastic sheets 14 , [l5| , 
growing tissues @, d, 0], and, possibly, graphene sheets 
[8| . Harnessing the elasticity of thin films artificially, 
however, remains both a theoretical and experimental 
challenge In recent work, Klein et al. 11 1 have 

designed thin films that swell isotropically, but inhomo- 
geneously. This creates a residual stress that is relieved 
by buckling out of the plane, leading to the possibility 
of controlling three-dimensional structure by patterning 
a two-dimensional surface. The problem of determining 
the optimal three-dimensional membrane shape from the 
underlying swelling pattern, however, remains a poorly 
understood but critical component of engineering func- 
tional thin membrane structures. 

For small thickness t, membranes exhibit large de- 
formations at relatively small cost compared to planar 
strain, and one generically expects that they will bend 
to eliminate most of that strain. Determining the re- 
sulting shape has been associated with the mathematical 
problem of isometric embedding P, 0, 0] _ how does one 
determine the buckled shape if one knows only how far 
apart points on the surface are supposed to be? When 
are the resulting shapes uniquely determined? When the 
buckling induces positive Gaussian curvature, it turns 
out that this is always possible and the buckled shape is 
unique [13j . For negatively-curved surfaces, on the other 
hand, the situation is murkier. Hilbert's theorem and its 
modern generalizations may forbid the vanishing of the 
strain globally even though the strain can always be ac- 
commodated by buckling locally [l]| . It is also known nu- 
merically that some isometric embeddings are not unique 
fisjl ■ However the precise conditions determining whether 
a strain-free embedding exists and is unique are not at 
all known in general. 
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FIG. 1: A typical disk with prescribed K = -1/R 2 (left) and 
an annular ribbon with three wrinkles (right), color-coded 
according to Gaussian curvature normalized to lie between —1 
and 1. The solid line on the ribbon is a curve of zero mean 
curvature. The 6-fold symmetry of the Gaussian curvature is 
apparent at the ribbon edge. 



In this letter, I study the physical manifestation of 
this geometrical frustration for thin, annular ribbons 
and disks that have been buckled by inhomogeneously 
swelling (for example, see figure[I]). The goal is to answer 
questions about the existence and uniqueness of buckled 
metrics with negative curvature in this restricted set of 
geometries, and determine how these properties physi- 
cally manifest themselves in buckled ribbons and disks. I 
will focus on the case of excess swelling on the outer edge 
of the ribbon resulting in a constant, negative Gaussian 
curvature. I find optimal ribbon shapes as the width in- 
creases, and identify a topological constraint that limits 
the availability of strain- free embeddings. These results 
are then applied to buckled disks with constant, negative 
curvature, where I present evidence for a critical radius 
above which the inhomogeneous strain can no longer be 
accommodated. Below the critical radius, the disk shape 
is uniquely determined by the metric. 
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GENERAL CONSIDERATIONS 

We consider a model thin film described by coordi- 
nates Xo in its flat, unswelled state. We assume either 
that there is a controlled growthprocess, such as cell divi- 
sion in a biological tissue [j, H, 0] , or the inhomogeneous 
swelling of a polymer gel [Uj , after which the film shape is 
given by X(xo) . The three-dimensional elastic energy can 
be described in terms of a Cauchy-Green strain tensor, 
A* ■ = 9X 1 /<9xq = djX 1 . We assume there is a unique 
energy minimum with respect to A 1 ■ . The question is to 
determine whether there is an X(xo) for which this A* 
exists. Consider the case of swelling an isotropic polymer 
gel. For the homogeneous swelling, A 1 - = fWj, were f2 
depends on the cross-link density as well as other ma- 
terial parameters [l2j]. When the cross- link density in 
inhomogeneous, as in Ref. [IH, I assume Q becomes of a 
function of the local density. For a thin film, dimensional 
reduction reduces the free energy to a two-dimensional 
membrane energy for the neutral surface of the film [l6| , 
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A (tr 7 ) 2 + 2 M tr ( 7 2 ) + k(H - c ) 2 , (1) 



where H is the mean curvature of the film, given by 
the average of the two principle curvatures, and 7 is 
the strain. In this letter, I will neglect a the possi- 
bility of spontaneous curvature, which would arise if 
£1 had a gradient across the film thickness, by setting 
co = 0. In general, the strain can be written in the 
form jij — — g^, where dl 2 = gijdx % dx^ is the two- 
dimensional metric of the surface and dl 2 = gijdx l dx^ 
is the prescribed two-dimensional metric. The trace is 
the usual one on curved surfaces, tr7 = g^jij, save that 
we use the prescribed metric for the contraction. Simi- 
larly, tr(7 2 ) = g Jijgi Jkl- For a membrane of thickness 
t, k/X oc t 2 . Therefore, the balance between bending 
and stretching energies is dominated by stretching and a 
reasonable approximation is to minimize the stretching 
energy completely, allowing the membrane to bend in 
whatever way necessary to do that. In general we should 
expect that the moduli vary with the swelling, but will 
neglect this complication in the following analysis. 

For isotropic swelling, the reference metric is dP = 
f2(r) [df 2 + f 2 d0 2 ] , where (f, 8) describe polar coordi- 
nates on the unbuckled membrane. We will find it use- 
ful, however, to choose a coordinate system in which 
the buckled membrane has a particularly simple mathe- 
matical representation. In Gaussian normal coordinates 
(r, 6), the metric will be dl 2 = dr 2 + p 2 (r,9)d9 2 , where 
p = r + lnf2(r) -1 ' 2 . These coordinates generalize the 
polar coordinates of a flat disk or ribbon: r measures 
the arc length along radial geodesies emanating from the 
center of the disk or inner ribbon boundary while 9 is the 
azimuthal angle around the disk or ribbon. These coordi- 
nates always exist in disk or ribbon geometries 13| . The 



boundaries are from r — — w/2 to r = tq + w/2 for 
ribbons and r = to r — w for a disk of radius w. The 
Gaussian curvature is given entirely in terms of the metric 
by K — —d 2 p/p. For positively curved disks, p increases 
slower than r; for negatively curved disks, it grows faster. 
For K = —1/R 2 , for example, p = i?sinh(r / R) which is 
linear near the origin and grows exponentially toward the 
edges. For the remainder of this letter, I will specialize 
to prescribed metrics with constant, negative curvature, 
though my methods will be applicable to a larger class 
of curvature distributions. 

Why can negatively-curved metrics not always be em- 
bedded? The classic no-go theorem was given by Hilbert, 
who showed that there is no surface of constant, nega- 
tive curvature that is both smooth and infinite area [1 31 ] - 
For axisymmetric shapes, this is apparently manifested 
physically by the divergence of the bending energy, which 
always occurs when d r p ~ 1 [9]. This is a compelling 
physical picture: the membrane is nearly strain-free un- 
til the curvature grows too large, after which the bending 
energy forces the film to have in-plane strain no matter 
how small the thickness. The situation, however, is more 
complex than it first appears. A surprising theorem of 
Nash demonstrated that a surface can be bent to have 
any metric but that the curvature may not be contin- 
uous (l9j |. In particular, this implies that the bending 
energy need not diverge in general [2(j. These theorems 
do not, unfortunately, lead to a general criteria one can 
use to determine when a negatively-curved metric can be 
strain-free. 

The key to our analysis is to examine the curvature 
more closely. The membrane curvature is described by 



the tensor hi 



didjT, where r(r,0) describes the 



embedding of the surface as a function of the coordi- 
nates, n is the unit normal to surface and the indices 
i = r,9. The mean curvature H = (l/2)(h rr + hgg/p 2 ) 
and (h rr hgg — h 2 e )/ p 2 — K is the Gaussian curvature, K. 
This formula for K, when combined with K — —d 2 p/p, 
yield a surprising result for isometric embeddings: the 
entire shape can be determined by integrating the evolu- 
tion equation 13l, |17| 



„ (n ■ d r dgv) - pd 2 p 
n ■ d„r 



(2) 



This equality, also known as Gauss' theorem egregium, is 
a second order evolution equation requiring initial condi- 
tions r(ro,0) and d r r(r ,8), since n = d r r x dgr/p. This 
initial data, however, is further constrained by compati- 
bility with the metric: (i) d r r 2 = 1, (ii) dgr 2 = p 2 and 
(iii) d r r ■ dgr — 0. This is sufficient to completely con- 
strain d r r up to an overall sign and, by equation ([2]), the 
rest of the surface as well. 

To determine the shape of wrinkled membranes, the 
evolution equation ^ can be recast as an evolution for 
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FIG. 2: HoR as a function of ro/R for different values of 
m. The thick line follows the mean curvature of the optimal 
ribbon of radius ro according to equation [5] up to m — 4. 
The transitions shown are at ro/R m 0.63 (m = 2 to 3) and 
r /R » 1.05 (ro = 3 to 4). 



the hij, yielding 



d r h r 



dgh r 



d r p 



d r hgg — h rr pd r p + dgh r g + 

pd 2 p = h 



d r p 



hgg — 



9gp, 



(3) 



hgghr 



Equations [3] are a form of the Gauss-Codazzi system [l3j |. 

If we use the last equation to solve for h rr we will have 
to contend with the fact that hgg may vanish. This is 
precisely what happens for axisymmetric surfaces, where 
hgg — > when d r p — > 1 implies a diverging mean curva- 
ture at one edge. This problem is exacerbated in disks, 
where <9 r /9 > 1 and no axisymmetric embedding exists. 
As it happens, however, hgg can vanish without H di- 
verging. Since 2H = [hf e + (hgg/ p) 2 - pd 2 p]/hgg, hnitc 
H requires that the components of the curvature tensor 
lie on a circle. In other words, 



h r g = \J p 2 H 2 + pd 2 pcosa, 
hgg/p = pH - y/ p 2 H 2 + pd 2 ps'ma, 
for some function a(r,6). This, in turn, implies 
h rr p = pH + y/ p 2 H 2 + pd 2 p sin a. 



(4) 



(5) 



This is akin to coordinate singularities in general rela- 
tivity; the problem is not with the surface but with the 
coordinate system on the surface. Notice the following 
interesting fact, however: we can choose any h rr and 
integrate equations (jSJ) to find the resulting metric with- 
out this difficulty. Apparently, the entire surface is de- 
termined by our choice of initial conditions and an arbi- 
trary function, h rr . The task will be to choose h rr such 
that the resulting p agrees with a prescribed metric p as 
closely as possible; as it turns out, this allows us to side 
step the issue of vanishing hgg. 




0.0 0.2 0.4 0.6 0.8 

rlR 



FIG. 3: Plot of normalized strain energy per unit radial 
length, £c(r) for a negatively curved disk with h rr increas- 
ing from zeroth order to eleventh order. The curves converge 
toward the rightmost curve with increasing order. 



ANNULAR RIBBONS 

I will first consider annular ribbons of width w and 
thickness t <C w. I will think of these ribbons as having 
been cut from a disk with metric p — Rsmh(r/R) from 
r = ro — w/2 to r = ro + w/2. This will greatly sim- 
plify our discussion of disks in the next section. Since 
the ribbons are thin, we solve the curvature evolution 
equations by expanding h rr — ^2 n h rr!Tl (5r / 'R) n and 
p = ^2 n p n {f/R) n - This will yield a unique, perturba- 
tive solution in powers of 5r = r — r for the shape. 

To make progress, we must make one further simpli- 
fication when the in-plane strain does not vanish. The 
strain contribution to the energy is controlled by A = 



(P 2 



2 )/p 2 ', after all, if A = then the actual metric 



equals the prescribed metric and 7^ — 0. Calculating the 
strain when A / can be quite difficult, however. If we 
assume A is small, we can expand all the strain compo- 
nents in powers of A and the dominant term in the strain 
energy will scale with A 2 . Clearly, minimizing A 2 is suf- 
ficient to find a strain-free configuration when one exists; 
I also expect that it provides an adequate approximation 
to the shape even when one does not exist. To obtain 
an explicit expression for the strain, I follow Audoly and 
Boudaoud [15| by setting all but one strain component 
to zero. I choose j rr = j r g = 0, leaving jgg = A. This 
approximation limits us, however, to regimes in which ei- 
ther in-plane strain or bending energy dominates the en- 
ergetics since it probably incorrectly captures how these 
two contributions to the energy would balance. In the 
case of buckled disks, we will only be concerned with the 
existence and shape of surfaces with A = so this should 
not be a problem. 
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As long asra>t the strain energy ensures that po — po 
and pi = pi - that is the metric agrees with the pre- 
scribed metric up to first order mw/R. The lowest order 
terms in the energy are then given by 
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where Y is the Young's modulus. Due to our approxima- 
tion for the strain, one should not play too close atten- 
tion to the numerical factors - I've included them here 
for completeness. The scaling of the compression and 
bending energies with respect to w/R allow us to iden- 
tify two regimes: w/R <C (t/R) 1 / 2 , for which the bending 
energy term dominates, and w/R^$> (i/i?) 1 / 2 , for which 
the compression term dominates. The former regime was 
studied in reference Q . Since the annular ribbons shapes 
are determined by minimizing their bending energy, they 
are flat (H(tq, 9) = 0) but circular strips. 

When w/R~^> (t/R) 1 / 2 , the dominance of the compres- 
sion strain results in pi = p2- After setting p2, the next 
term in the compression energy scales as (w/R) 7 . More 
generally, if we choose p„ = p n , the dominant term in 
the compression energy scales as (w/R) 3+2n . From this, 
we identify a cascade of regimes (t/R) <C (t/R) 1 / 2 <C 
(t/R) 1 / 4 <?;••• that set the relative balance of compres- 
sion energy and bending energy order by order. When 
w/R ^> (t/R) 1 / 4 , for example, the compression energy 
requires us to choose p$ = p% . We then obtain an approx- 
imation to the ribbon shape by minimizing the compres- 
sion strain up to some power of w/R, then set the remain- 
ing terms by minimizing the bending energy. This basic 
picture will not be changed by a more careful evaluation 
of the strain energy, and if the ribbons were not closed 
annuli the problem of optimal shape would be solved. 
In the limit of t — > 0, in particular, we conclude that 
we should systematically set the coefficients p n = p n in 
the ground state. This analysis implies that we are able 
to choose any H and a arbitrarily for our initial con- 
ditions, from which the entire expansion in powers of r 
is determined. This analysis conforms with the notion 
that there is no local obstruction to finding strain-free 
configurations. 

For annular ribbons, however, there is an obstruction 
to our ability to minimize the compression energy term to 
third order and above: the ribbon must be periodic and 
not twisted. This constraint is encoded by the winding 
numbers of the three Euler angles of the coordinate frame 
defined by d r r, dgr and n. The unit tangent vector of a 
curve of constant r can be written as t = cos 9 sin cpx — 
cos 9 cos <py — sin^z. The remaining Euler angle, ip de- 
fines angle made by the surface normal n around the axis 
defined by t: n — (— sin?/ cos + sin0cos , (/'sin</))a; + 
(— sin?/! sin — sin 9 cos tjj cos 9)y + cos 6* cos -02:. Using 



these definitions, we have sin ip9' — cos 9 cos ipft = —p\ 
from the metric, where the prime indicates differentia- 
tion with respect to the azimuthal coordinate. These 
angles are related to the curvatures by 



tp = tan# [pi cos"0 — smiphee 
9' = — cosiphee/ Po — sinippi 
, cos^pi - sin?/) hee/po 



/Po] 



(7) 



cos 9 

The constraint states simply that A(j> = 0(2-7r) — 0(0) = 
2ir, Atp = 0, and A9 = 0, or some other combination 
where one Euler angle gives 2tt and the others give zero. 

We look for a solution of these equations of the form 
H = _ffosin(m#), a = m9 and p2 = p2 in equation l|4|). 
choosing the number of wrinkles m, and Hq, to satisfy 
the winding number constraints. The numerically deter- 
mined relationship is shown in figure O Apparently, Hq 
is determined completely by m and p\ , and it is no longer 
possible to set the curvature to zero. Moreover, this re- 
lationship must be true not only for ribbons but for all 
curves of constant r on a membrane. Notice that Hq = 
along r = tq whenever m = p\ + 1. Figure [T] displays a 
ribbon with m = 3 and Hq — at r = r$. 

With these results in hand, we consider again the case 
that w/R^> (t/R) 1 / 4 in which we are required to choose 
a minimal value of p^. The evolution equation for p de- 
termines p3 in terms of the curvatures at some r = tq. 
Assuming HqR is small, we expand the compression en- 
ergy to quadratic order in the mean curvature, yielding 
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This expression can be minimized subject to the topo- 
logical constraint. The results of this minimization are 
displayed as a thick solid line in figure [2] in terms of 
the mean curvature. In particular, I find that m = 2 
is preferred for small radii, m — 3 is energetically pre- 
ferred when r w 0.63i?, and m — 4 when r ~ 1.05i?. 
This trajectory differs from the optimal ribbon in the 
regime (t/R) 1/2 < w/R < (t/R) 1/4 , when bending en- 
ergy dominates over compression strain. In that case, the 
number of wrinkles follows the trajectory in figure [5] that 
minimizes H 2 . This fundamental incompatibility arises 
from the topologically determined relationship between 
Hq and the metric - if we were to cut the ribbon to relax 
this constraint, a ground state with Hq = and ps = 
would yield a non-integer m according to equation [5J 



BUCKLED DISKS 

We now turn our attention to disks. Since every curve 
of constant r on a disk satisfies the same topological con- 
straint that ribbons do, we can recover the limit of a disk 
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FIG. 4: HR for the disk of figure EE] for r = R,3R/2, and 27?. 

by taking r$ — > in figure [2 This implies that p\ — > 1 . 
Notice that when to = 2, i? (0) = and the disk is smooth 
at the center. On the other hand, if to > 2 on the disk 
as ro — > 0, the surface will not be smooth at r = since 
H = Ho(0)sm(m8) - the result is a singularity at the 
disk center akin to that of a cone. 

We again apply an expansion in powers of r/R. We 
will assume that the t is small enough that the com- 
pression strain dominates to all orders of r/R. This 
also ensures that the resulting shape will be indepen- 
dent of the approximation used to evaluate the strain 
since we need only minimize A 2 . Immediately, we see 
that p = Y,nPn(Sr/R) n agrees with p = J2 n P~n($r/R) n 
up to second order when h rr is given by equation ([3]) 
evaluated on r — tq. Using this expansion, I find that 
a = 29 as ro — > as expected from our analysis of figure 
[2j Negatively-curved disks are saddle-shaped near their 
centers regardless of the magnitude of Gaussian curva- 
ture. This result is anticipated by a stability analysis of 
growing tissues Q and the fact that saddles are canon- 
ical examples of negatively-curved surfaces. Continuing 
in this manner, 

h rr = [l + (l/2)5r 2 /R 2 + (l/6)5r 4 /R i ]R- 1 sm(28) 
-{l/24)5r 4 /R 4 sm(69) (9) 

up to fourth order, implying that p agrees with the pre- 
scribed metric up to seventh order. The resulting shape 
is shown in figure [1] It is worth noting that for small 
but finite t/R, we might expect deviations from this 
shape due to small, but nonzero, contributions from the 
bending energy. In this case, the resulting deviations of 
the Gaussian curvature are fourfold symmetric. In fact, 
K oc sin 2 (29), which is a consequence of its quadratic de- 
pendence on the hij and the membrane's saddle shape. 
This four-fold symmetric deviation is apparent in exper- 
iments on non-Euclidean disks (ll| and is universal for 
negatively-curved disks with axisymmetric p. 

This perturbative construction can be continued to 
higher orders, forcing the introduction of higher oscil- 
lations in 9. The data in figure [3] is shown up to eleventh 



order in r/R, after which the deviation of the metric ap- 
pears to converge to a universal curve. I have confirmed 
this convergence up to fifteenth order in r/R. It is inter- 
esting that the perturbative expansion is quite good, as 
measured by A 2 , even with relatively few terms, as long 
as r < R. From figure [31 above a radius, r ss 1.2R, A 2 
does not decrease with increasing order of r/R though it 
does below this radius. 

It is enlightening to consider these results in light 
of several rigorous mathematical results on negatively 
curved surfaces. Since a surface can be deformed to agree 
with any prescribed metric in a finite region 17[ we con- 
clude that the unique expansion in equation §§§ must 
have a finite (nonzero) radius of convergence around the 
origin. Hilbert's theorem, however, ensures that this ex- 
pansion fails to converge at large enough r even at t — > 
since, if it did, the surface would have infinite area and 
constant negative curvature. A more recent mathemat- 
ical result does provide an existence proof (but not a 
construction) that a negatively-curved disk can be em- 
bedded up to a critical radius r c lj, LL8J . According to 
this result, above this radius \K\ must either decrease 
faster than 1/r 2 or develop zeros no matter how h rr is 
adjusted. From figure [3j I find r c « 1.2R as the criti- 
cal radius for constant curvature. These theorems, taken 
together, imply the disk becomes non-analytic at r = r c 
as t — » 0. This follows immediately from the fact that 
p 2 — p 2 — for r < r c and is nonzero for r > r c . 

Notice that figure [4] indicates that H = Ho(r)s'm(28) 
on a disk to excellent approximation. We glean fur- 
ther insights by considering the topological constraint 
on curves of constant r on the disk. As it turns out, 
J d9Ho(r, 9) 2 /(2ir) on the disk agrees quite well with the 
value predicted by the to = 2 branch of figure[2] This pre- 
diction works until the m — 2 branch ends at r ~ 1.2R, 
beyond which the constraint cannot be satisfied by solu- 
tions of the form H — Hq sin(20). This gives an indepen- 
dent corroboration of the value of critical radius found by 
perturbative expansion in powers of r/R. Unlike the case 
of axisymmetric surfaces, the mean curvature does not 
diverge at this critical radius. This roughly matches our 
intuition from Nash's theorem that the curvature need 
only be discontinuous at worst, not divergent. It is in- 
teresting to speculate that the curvature would become 
discontinuous precisely at r = r c as t — ► if we pursued 
an isometric embedding beyond the critical radius. 



CONCLUSION 

What happens above the critical radius? A hint comes 
from the fact that better agreement between disk shape 
and prescribed metric can be made if we allow additional 
subwrinklcs. In figure El I compare the constant cur- 
vature saddle to a disk on which additional subwrinkles 
begin abruptly at r/R = 4/3 with m — 10. This was cho- 
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FIG. 5: p(r) for the saddle (up to order r 4 ) with and without 
wrinkled edges, compared to the prescribed metric p(r) = 
Rsinh(r / R) (dashed). For comparison, p(r) — r (dotted), a 
flat disk, is shown. 



should be directly applicable to more general metrics. 
Future work will focus on the defining the precise con- 
ditions under which strain-free shapes exist, understand- 
ing which metrics uniquely specify a shape, and mak- 
ing a more direct connection between experimentally- 
controlled parameters such as the cross-link density and 
the buckled film shape. Finally, the predictions made by 
this theory can be tested on thin films with prescribed 
inhomogeneous swelling patterns. 

I am grateful to G. Grason and R.D. Kamien for useful 
discussions and a reading of the manuscript. I thank the 
Aspen Center for Physics, where a portion of this work 
was done. 



sen to illustrate the significant benefit in strain achieved 
by adding subwrinkles. Subwrinkling is often observed 
in negatively-curved leaves @, 0] , strips 14 1 , flowers Q , 



and disks It is tempting to speculate that the zero 



compression strain trajectory of figure [2] also indicates 
the preferred number of wrinkles as a function of radius 
on a disk. Unfortunately, this does not account for the 
energy cost of changing the number of wrinkles and so it 
is unlikely that disks will exhibit all the m as they do for 
ribbons. Further insights will be needed. 

In summary, I have introduced new methods to com- 
pute optimal shapes of annular ribbons and thin disks. 
This theory identifies a cascade of shape regimes as a 
function of ribbon width and makes definite predictions 
for the resulting optimal shapes. In particular, I con- 
sider the new regimes (t/R) 1/2 < w/R <C (t/i?) 1/4 and 
(t/R) 1 / 4 <C (w/R) <C (t/R) 1/e . The novel feature in 
these regimes are the interplay between topology and 
curvature, which constrains the number of wrinkles and 
mean curvature. These topological constraints also apply 
to closed curves on thin disks, and can be used as good 
approximations to the shape of negatively curved disks. 

One obvious limitation of the theory is the computa- 
tion of the strain energy, which sets two strain compo- 
nents to zero. This is most assuredly incorrect. In the 
present work, I have limited consideration to regimes in 
which there is a clear dominance of either the compres- 
sion or the bending energy order by order in perturba- 
tion theory. Therefore, minimization of the compression 
strain amounts to setting the appropriate power of the 
metric equal to its prescribed value. This assumption 
breaks down beyond the critical radius r c on a disk, as 
well as in the cross-over regions between two regimes of 
an annular ribbon - for example, when w/R ~ (t/R) 1 / 4 
- where compression and bending energies are roughly 
balanced. The ribbon shape in these cross-over regions 
remains elusive. 

Though this letter is limited to disks and ribbons with 
constant negative curvature, the methods of analysis 
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